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Paring correlations in weakly bound halo nuclei 6 He and 11 Li are studied by using a threc- 
■ body model with a density-dependent contact interaction. We investigate the spatial 

' structure of two-neutron wave function in a Borromean nucleus n Li. The behavior of 

the neutron pair at different densities is simulated by calculating the two-neutron wave 
function at several distances between the core nucleus 9 Li and the center of mass of the 
two neutrons. With this representation, a strong concentration of the neutron pair on the 
nuclear surface is quantitatively established for neutron-rich nuclei. Dipole excitations 
^ ' in 6 He and 11 Li are also studied within the same three-body model and compared with 

C _ ) experimental data. The small open angles between the two neutrons from the core are 

extracted empirically by the B(E1) sum rule together with the rms mass radii, indicating 
C^") ' the strong di-neutron correlation in the halo nuclei. 
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Pairing correlations play a crucial role in many Fcrmion systems, such as liquid 3 He, 
atomic nuclei, and ultracold atomic gases. When an attractive interaction between 
two Fermions is weak, the pairing correlations can be understood in terms of the 



well-known BCS mechanism, that shows a strong correlation in the momentum 
space. If the interaction is sufficiently strong, on the other hand, one expects that 
two Fermions form a Bosonic bound state and condense in the ground state of many- 
body system. 1 The transition from the BCS- type pairing correlation to the Bosc- 
Einstein condensation (BEC) takes place continuously as a function of the strength 
of attractive interaction. This feature is referred to as the BCS-BEC crossover. 

It has been feasible by now to study the structure of nuclei on the edge of 
neutron drip line. Such nuclei are characterized by a dilute neutron density around 
the nuclear surface so that one can investigate the pairing correlations at several 
densities, 2 ranging from the normal density in the center of nucleus to a diluted 
density at the surface. The pairing correlations are predicted to be strong at the 
surface, but rather weak at the normal density and also far outside of the core. Thus, 
the weakly bound nuclei will provide an ideal environment to study the dynamics 
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of pairing correlations in relation with the BCS-BEC crossover phenomenon. 

In this talk, we discuss the manifestation of the BCS-BEC crossover phenomenon 
in finite neutron-rich nuclei. We particularly study the ground state wave function 
of a two-neutron halo nuclei, 6 He and 11 Li. These nuclei are known to be well 
described as a three-body system consisting of two valence neutrons and the core 
nucleus( 4 He or 9 Li). 3 ~ 5 A strong di-neutron correlation as a consequence of pairing 
interaction between the valence neutrons has been shown theoretically in 11 Li. 5 ' 6 
We take this nucleus to study BCS-BEC crossover features in connection to the 
strong two-neutron correlation, which has recently been observed experimentally in 
low-lying dipole strength in 11 Li. 7 



2. Three-body model and di-neutron correlation 

In order to study the pair wave function in n Li, we solve the following three-body 
Hamiltonian, 5 ' 6 

H = h nC (l) + h nC {2) + V nn + ^-^ 1 (1) 

J\ C TTL 

where m and A c are the nucleon mass and the mass number of the inert core 
nucleus, respectively. h n c is the single-particle Hamiltonian for a valence neutron 
interacting with the core. We use a Woods-Saxon potential for the interaction in 
h n c- The diagonal component of the recoil kinetic energy of the core nucleus is 
included in h n c, whereas the off-diagonal part is taken into account in the last term 
in the Hamiltonian (TT|). The interaction between the valence neutrons V nn is taken 
as a delta interaction whose strength depends on the density of the core nucleus. 
Assuming that the core density is described by a Fermi function, it reads 

V nn (n,r 2 ) = S(n - r 2 ) (v + — Vp ) , (2) 

V 1 + exp[(i? - R p )/a p \J 

where R = + 7^*2 ) / 2 1 . We use the same value for the parameters as in Refs. 5,6 

The two-particle wave function ^(fi, r-i) is obtained by diagonalizing the three- 
body Hamiltonian (fTJ) with a large model space which is consistent with the nn 
interaction, V nn . To this end, we expand the wave function \t(fi, r%) with the eigen- 
function 4> n ijj z (ri) of the single-particle Hamiltonian h n c- In the expansion, we 
explicitly exclude those states which are occupied by the core nucleus. 

The ground state wave function is obtained as the state with the total angular 
momentum J = J z — 0. We transform it to the coordinate system with the relative 
and center of mass (cm) motions for the valence neutrons, r = f 1 — f 2 and R = 
(fi + r*2)/2. 8,9 The wave function is first decomposed into the total spin 5=0 and 
S=l components. Then, the coordinate transformation is performed for the S=0 
component, which is relevant to the pairing correlation: 

* S=O (fi,30 - £ h(r, R) [Y L (f)Y L (k)}^ \ Xs =o), (3) 

L 
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where \\s=o) is the spin wave function. The two-particle wave function is plotted 
for n Li in two different coordinates in Fig. [TJ Fig. [lja) is plotted as a function of 
r = r\ — T2 and the angle between the valence neutrons, while the radial coordinates 
r and R are adopted for Fig. [Tfc). The S — component is only taken for Fig. [TJb). 
One observes two peaks in both the figures. The peak at smaller angle 0\2 in Fig. 
QJa) is referred to as "di-neutron configuration" , while the large angle is called 
"cigar-like configuration" . One can sec the di-ncutron configuration has a long tail 
as a typical feature of halo wave function. Figure Q] (c) shows the square of two- 
particle wave function for the L — component. One can clearly recognize the two 
peaked structure in the plot, corresponding to the di-neutron and the cigar-like 
configurations. 

The L = wave functions of 11 Li for different values of R are plotted in Fig. 
El Since we consider the density-dependent contact interaction, Eq. ([2]), this is ef- 
fectively equivalent to probing the wave function at different densities. At R = 0.5 
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Fig. 1. A two dimensional plot of the square of the ground state two-particle wave function 
for xl Li; (a) the total density, (b) the S = component as a function of the radial coordinate 
r = ri = r'2 and the angle between the valence neutrons 812- (c) f 2 R 2 |/^ = o(r, R)\ 2 in Eq. JSJ, as 
a function of the relative distance r and the center of mass coordinate R for the valence neutrons 
as denoted in the upper right inset. 
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Fig. 2. The ground state two-particle wave functions, r 2 R 2 \fj J= o(r, R)\ 2 of 11 hi as a function of 
the relative distance between the neutrons, r, at several center of mass distances R as indicated 
in the inset. Notice the different scales on the ordinate in the various panels. 

fm, where the density is close to the normal density po, the two particle wave func- 
tion is spatially extended and oscillates inside the nuclear interior. This oscillatory 
behavior is typical for a Cooper pair wave function in the BCS approximation, and 
has in fact been found in nuclear and neutron matters at normal density po- 10 ' 11 As 
R increases, the density p decreases. The two-particle wave function then gradually 
deviates from the BCS-like behavior. At R = 3 fm, the oscillatory behavior almost 
disappears and the wave function is largely concentrated inside the first node at r ~ 
4.5 fm. The wave function is compact in shape, indicating the strong di-neutron cor- 
relation, typical for BEC where many such pairs are present. At R larger than 3 
fm, the squared wave function has essentially only one node, and the width of the 
peak gradually increases as a function of R. This behavior is qualitatively similar 
to the pair wave function in infinite matter. 10 We have confirmed using the same 
three-body model that this scenario also holds for another Borromean nucleus 6 He 
as well as for non-Borromean neutron-rich nuclei 16 C and 24 0. 

The transition from the BCS-type pairing to the BEC-type di-neutron correla- 
tion can also clearly be seen in the root mean square (rms) distance of the two neu- 
tron system. We plot this quantity in Fig. 3(a) as a function of R. In order to com- 
pare it with the rms distance in nuclear matter, we relate the cm distance R with the 
density p using the Fermi-type functional form p(R)/po = [l + exp((iZ — R p )/a p )] , 
as used in the nn interaction in Eq. ([2|). Fig. 3(b) shows the rms distance as a func- 
tion of density p thus obtained. The rms distance shows a distinct minimum at 
p ~ 0.4po (R ~ 3.2 fm). This indicates that the strong di-neutron correlation grows 
in 11 Li around this density. Notice that the probability to find the two-neutron pair 
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Fig. 3. The root mean square distance r rma for the neutron pair. 

is maximal around this region (see Fig. [1] (c)). The behavior of rms distance as a 
function of density p qualitatively well agrees with that in infinite matter (see Fig. 
3 in Ref. 10 ), although the absolute value of the rms distance is much smaller in the 
finite nucleus. 



3. Dipole excitations and correlation angles 



The rms distance J (rj?_ 2r ,) has an intimate relation to the B(E1) strength as, 3,4,12 

B{El) = l(^X (rl_ 2n ). (4) 



n \ A , 

This relation is obtained with closure, which includes unphysical Pauli forbidden 
transitions to the states with negative excitation energies. Although the effect of 
Pauli forbidden transitions is not large, it leads to a non-negligible correction. In 
Ref., 12,13 it has been proposed to estimate the experimental value for (r 2 _ 2jl ) using 
the relation, 

, 2 \ ._ B(El;E < E max ) cxp 2 

Vc-2n/exp— „, j-, „ ^ \' c _2n/cal- V J 7 

D\£j\.,hj <. £/maxJcal 

The dipole strength distributions for the 6 He and 11 Li nuclei obtained with the three 
model are shown in Fig. 0] Also shown by the solid curves are the B(E1) distributions 
smeared with the Lorenzian function with the width of T = 0.2 McV. For the 6 He 
nucleus, we obtain the total B(E1) strength of 0.660 e 2 fm 2 up to E < 5 MeV and 
1.053 e 2 fm 2 up to E < 10 MeV. These are in good agreement with the experimental 
values, B(E1; E < 5 MeV)=0.59 ± 0.12 e 2 fm 2 and B(E1; E < 10 MeV) = 1.2 ± 0.2 
e 2 fm 2 . 15 For the l: Li nucleus, we obtain the total B(E1) strength of 1.405 e 2 fm 2 up 
to -B re i = E — S2n < 3 MeV, which is compared to the experimental value, B(E1; 
E Te i < 3 MeV)=1.42 ± 0.18 e 2 fm 2 . 7 Again, the experimental data is well reproduced 
within the present model. From the calculated values for (? ,2 _2n)cab that is, 13.2 and 
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26.3 fm 2 for 6 He and n Li, respectively, we thus obtain J{r 2 c _ 2n ) exp = 3.878± 0.324 
fm and 5.15 ± 0.327 fm for 6 He and 11 Li, respectively. Notice that the value for the 
6 He nucleus is somewhat larger than the one estimated in Ref., 15 that is, 3.36 ± 
0.39 fm. 

We next evaluate the Coulomb breakup cross sections based on the relativistic 
Coulomb excitation theory. 1 These are obtained by multiplying the virtual photon 
number N-ei(E) to the B(E1) distribution shown in Fig. 0] The solid line in Figs. 
[5] (a) and (b) shows the Coulomb breakup cross sections thus obtained for 6 He+Pb 
reaction at 240 MeV/nucleon 15 and 11 Li+Pb reaction at 70 MeV/nucleon, 7 respec- 
tively. In order to facilitate the comparison with the experimental data, we smear 
the discretized cross sections with the Lorenzian function with an energy depen- 




Fig. 5. Coulomb breakup cross sections (a) for 6 He+Pb at 240 MeV /nucleon, and (b) for 
11 Li+Pb at 70 MeV /nucleon. The solid line is the result of the full three-body calculations, 
while the dashed and the dotted lines are obtained by treating the recoil term approximately (see 
Ref. 13 for details). These results are smeared with an energy dependent width of T = 0.15 ■ \ZE IC \ 
MeV. The experimental data are taken from Ref. 15 for 6 He and from Ref. 7 for 11 Li. 
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dent width, T — a ■ y/E^h We take a = O.lSMeV 1 / 2 and 0.25 MeV 1 / 2 for 6 He and 
11 Li, respectively. We see that the experimental breakup cross sections are repro- 
duced remarkably well within the present three-body model, especially for the n Li 
nucleus. 

Let us now discuss the geometry of the 6 He and n Li nuclei. Using the experi- 
mental value for (r^,_ 2n ) obtained from the B(E1) distribution, one can extract the 
mean opening angle between the valence neutrons once an additional information is 
available. The mean opening angle can be extracted directly when the rms distance 



between the valence neutrons, 
matter radius and 



is available. This quantity is related to the 



-In i 



in the three-body model, 3 ' 6 ' 12 



-4, 



2A C 



( r m) — ~T ( r m)A c + — PT 



1 



vmi^-c • v c-lnl ■ 

where A c = A — 2 is the mass number of the core nucleus. The matter radii (r 2 ^) 
can be estimated from interaction cross sections. Employing the Glauber theory 
in the optical limit, Tanihata et al. have obtained */ (r 2 n ) = 1.57 ± 0.04, 2.48 ± 
0.03, 2.32 ± 0.02, and 3.12 ± 0.16 fm for 4 He, 6 He, 9 Li, and n Li, respectively. 16 
Using these values, we obtain the rms neutron-neutron distance of \J (r 2 ra ) = 3.75 
± 0.93 and 5.50 ± 2.24 fm for 6 He and 11 Li, respectively. Combining these values 
with the rms corc-di-ncutron distance, yj (r 2 _ 2 n): obtained with Eq. |5]), we obtain 

the mean opening angle of (9 nn ) = 51.561^2 2 aim 56.2^2i'3 degrees for 6 He and 
11 Li, respectively. These values are close to the result of the three-body model 
calculation, (6 nn ) =66.33 and 65.29 degree for 6 He and 11 Li, respectively, 5 although 
the experimental values are somewhat smaller. An alternative way to extract the 
value a/ (r 2 „) was reported by the three-body correlation study in the dissociation 
of two neutrons in halo nuclei. 17,18 The two neutron correlation function provides 
the experimental values for J (r 2 n ) to be 5.9 ± 1.2 and 6.6 ± 1.5 fm for 6 He, 11 Li, 

respectively. 17 When one adopts the presently obtained value for y (r 2 _ Q „) with 

Eq. ([5]) instead of those in Refs., 7,15 one obtains 



(6) 



>=74.5I^ and 65.2 + 13 . 



2n/ 

+ 11 ' 4 for 



6 He and 11 Li, respectively. Notice that these values are in good agreement with the 
results of the three-body calculation 5 as is seen in Table 1. 



Table 1 . The geometry of the 6 He and 1 1 Li nuclei extracted from various experimental 
data. The mean opening angles calculated by the three-body model are also given in 
the last line for each nucleus in the table. 



nucleus 


yj(rl_ 2n ) (fm) 


V(r*„) (fm) 


method 


(Onn) (deg.) 


6 He 


3.88±0.32 


3.75 ± 0.93 
5.9± 1.2 


(matter radii) 
(neutron correlations) 


51 56 +11 -' 2i 

Oi.OD _ 12 .37 

74 5 + 11 ' 2 

(4.0 _ 131 

66.33 s 


"Li 


5.15±0.33 


5.50 ± 2.24 
6.6± 1.5 


(matter radii) 
(neutron correlations) 


56.2 +J<-» 
65 2 + 11 ' 4 

DO.Z _ 13 o 

65.29 s 
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4. Summary 

We studied the two-neutron (2n) wave function in the Borromean nuclei 6 He and 
11 Li by using the three-body model with the density-dependent pairing force. We 
explored the spatial distributions of 2n wave function as a function of the cm dis- 
tance R from the core nucleus and found that the structure of the 2n wave function 
alters drastically as R is varied. We also showed that the relative distance between 
the two neutrons scales consistently to that in the infinite matter as a function of 
density. These features are in close analogue to the characteristics of the BCS-BEC 
crossover phenomenon found in the infinite nuclear and neutron matters. We have 
used the same three-body model to analyze the B(E1) distribution as well as the 
Coulomb breakup cross section of the 6 He and 11 Li nuclei. We have shown that the 
strong concentration of the B(E1) strength near the continuum threshold can be 
well reproduced with the present model for both the nuclei. Using the calculated 
B(E1) strength, we extracted the experimental value for the rms distance between 
the core and di-neutron, which was then converted to the mean opening angle of the 
two valence neutrons. We have found that the mean opening angles thus obtained 
are in good agreement of the results of the three-body model calculation. 
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